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Abstract The discovery of dark energy by the first generation of high-redshift supernova 
surveys has generated enormous interest beyond cosmology and has dramatic imphcations for 
fundamental physics. Distance measurements using supernova explosions are the most direct 
probes of the expansion history of the Universe, making them extremely useful tools to study 
the cosmic fabric and the properties of gravity at the largest scales. The past decade has seen 
the confirmation of the original results. Type la supernovae are among the leading techniques 
to obtain high-precision measurements of the dark energy equation of state parameter, and in 
the near future, its time dependence. The success of these efforts depends on our ability to 
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understand a large number of effects, mostly of astrophysical nature, influencing the observed 
flux at Earth. The frontier now lies in understanding if the observed phenomenon is due to 
vacuum energy, albeit its unnatural density, or some exotic new physics. Future surveys will 
address the systematic effects with improved calibration procedures and provide thousands of 
supernovae for detailed studies. 
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1 Introduction 

Particle physics and cosmology are tightly linked and ideally complement each 
other. The physics of the early universe is beyond the reach of even the most pow- 
erful accelerators and needs to be deduced from the embers left by the Big Bang. 
The richness of the cosmic microwave background radiation as a source of infor- 
mation of the early universe has been described many times and has to be counted 
amongst the most successful scientific endeavors of the past two decades [THS]. 
Observational cosmology can now address particle physics questions, e.g. the 
number of neutrino species and the sum of their masses, beyond the Standard 
Model particles as dark matter candidates, and even braneworld scenarios in- 
voking extra spatial dimensions. Theoretical high-energy physics has started to 
address questions concerning the uniqueness of the observable universe and String 
Theory inspired models to explain the current state of the cosmos [HIS]. Super- 
novae, the violent destruction of an entire star, offer means to obtain accurate 
information about the expanding universe, and are therefore a precious instru- 
ment in the cosmology tool box. These luminous explosive events are important 
for modern cosmology since they can be detected even in very remote galaxies. 
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Furthermore, some subtypes, Type la supernovae (SNe la) in particular, can be 

calibrated fairly reliably to provide accurate distances that can be used to map 

the expansion history of the universe. 

The current expansion rate, the Hubble constant Hq, was determined with su- 
pernova distances for several decades [BHH]. Two different methods have been em- 
ployed using supernovae. Assuming a unique luminosity for SNe la, often referred 
to as the "standard candle" method, one can directly use the observed bright- 
ness to infer the distance, an appealing technique due to its simplicity. Various 
(astrophysical) effects influence the distance measurement, but through suitable 
calibration and corrections this remains very successful. The other method is to 
measure the physical expansion of a Type II supernova through the radial veloc- 
ity and its brightness increase reflecting the growth in surface area. This method 
has been applied to supernovae stemming from massive stars and with extended 
envelopes. In the local universe, this measurement is still possible, but due to the 
limited luminosity of the Type II supernovae becomes very difficult with current 
observing facilities at larger distances [9]. 

Dark Energy, detected using SNe la through the accelerated expansion of the 
universe, is amongst the most notorious recent additions to our understanding 
of the cosmic composition (e.g. [TUHHj). After the first indications based on 
only small samples of supernovae [T5l[16] the past decade has seen a significant 
amount of telescope time invested in further refining our understanding of SNe la 
and their ability as distance indicators. The cosmological results have been con- 
firmed and the uncertainties have decreased dramatically due to better controlled 
observational techniques, improved calibrations, larger samples and an extended 
redshift range. The current studies are limited by the systematic uncertainties re- 
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lated to the supernova explosions and the hght propagation through the universe 

rather than the size of the statistical sample. Overall, there are now several hun- 
dred Type la supernovae available for cosmic distance measurements and there 
are prospects to increase these numbers at least ten fold within the next decade. 
The field has moved from the discovery of dark energy to its first characteriza- 
tion through the measurement of the equation of state parameter, w = p/pc^, i.e., 
the ratio between pressure and energy density when treated as a fluid. Future 
projects will address potential time-dependence of this parameter, w{t). Note 
that from now on, we will use natural units, i.e., c = 1, unless stated otherwise. 

Although there is no observation today that is incompatible with Einstein's 
cosmological constant. A, being responsible for the accelerated expansion of the 
universe, only few high-energy physics theorists are ready to quit searching for 
alternative models. 

A cosmological constant {w = —1) can be naturally associated with the vacuum 
zero-point energy predicted by quantum field theory. However, the corresponding 
density, p^ac ~ (10~^eV)^, does not match any particle physics scale, and is more 
than 120 order of magnitudes lower than the "natural" scale set by the Planck 
mass, Mp ~ (lO^^eV)^or about 60 orders of magnitude below the supersymmetry 
breaking scale, M^jj^y ~ (lO^^eV)^. Furthermore, it may be regarded as a great 
coincidence that p^ac and pM, are so close in spite of pM oc a~^, i.e., the density 
of non-relativistic matter gets diluted proportionally to the growing volume of 
the expanding universe, while the vacuum energy density remains constant. 

Some highly regarded physicists have invoked the Anthropic principle to elude 
these difficulties arguing that we are witnessing the result of a selection process: 
galaxies, stars, planets and intelligent beings would only have a chance to come 
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to being in "un-natural" universes like ours [171118) . Another common view on 

this problem is that a yet unknown mechanism sets the vacuum energy to zero 

and the present day acceleration is instead caused by new physics, e.g., a light 

scalar field or an effective change of gravity at the largest scales in a higher 

dimension universe. While the Anthropic hypothesis cannot be further tested 

experimentally (with a possible exception described in [19j), we will not address 

it further here. Alternatives to A, however, would imply that the dark energy 

density would change in time (and space), something that can be tested with 

SNIa, especially when combined with other probes. In particular, measurements 

of the expansion history and growth of structure in a universe with dynamically 

evolving dark energy would lead to 'w{t) ^ —1, i.e., distinguishable from A. 

A new generation of multi-probe projects are being planned, and their relative 
merits have been assessed by the ESA-ESO Working Group on Fundamental 
Cosmology [20] and the Dark Energy Task Force report [21J. These forthcoming 
measurements are expected to provide the first accurate measurements of the 
time dependence of w. 

The use of supernovae for cosmology has been reviewed many times. The 
original cosmological application of supernovae was for the Hubble constant [71|8l 
I22j . The connection of supernovae with Dark Energy has been described in many 
reviews as well |10 U 12 1[23U25j . Using supernovae to map the cosmic distance scale 
requires a sufficient understanding of the explosions and the astrophysical effects, 
which affect light propagation. These are very important considerations for the 
derivation of accurate distances. 
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1.1 Supernova classification 

Supernovae have been recognised as an own class of astrophysical objects by 
Baade and Zwicky following earlier work by Lundmark [27], who noted that 
there is a class of eruptive variable objects, which are about 10 magnitudes (i.e. 
10000 times) more luminous than regular (and more frequent) novae. An early 
classfication scheme for supernovae was introduced by Minkowski [281129] . 

Supernova classification is based on their spectral appearance near maximum 
light [30j. The Type II supernovae show hydrogen in their maximum-light spec- 
trum, while Type I supernovae lack hydrogen. Type lb SNe display prominent 
helium lines, while the Type Ic display neither hydrogen nor helium. Type la su- 
pernovae spectra are dominated by lines from higher-mass elements, like calcium, 
sulphur, silicon and iron, but lack hydrogen and helium - the most abundant el- 
ements in the universe. 

The past years have seen some interesting modifications to the classification 
after several decades of relative few new additions [29l|31]. Several examples of 
objects, which do not fit any of the established classes have been added. The Type 
la class is overall very uniform in its appearance, but shows some spread in its 
luminosity and spectral evolution. Both overluminous and underluminous SNe la 
have been identified (e.g. P^). In particular, a range of expansion velocities - 
as derived from the absorption trough of several lines (e.g. |32fl35j ) has been 
established. These peculiar objects typically show interesting, and sometimes 
subtle, deviations in their spectral appearance, which may allow to refine their 
use as distance indicators (see Sec. [5]). Thanks to the large number of SN surveys 
conducted over the past decade, a small number of truly different objects has 
been discovered. Among them are objects, which appear very close to SNe la. 
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but displayed a strong hydrogen emission (SN 2002ic; [36]), extremely luminous 

objects with rather slow expansion velocities (SN 2003fg; [37J, SN 2006gz; [38] . 
SN 2009dc; [39|), and extremely faint objects (SN 2000cx; ^40], SN 2002cx; 
2005hk; [32], 2008ha; [33], but see [Hj for a differing view on the classification of 
this object). In many cases, the peculiar nature of these objects only became clear 
through a detailed analysis of the pre-maximum color and spectral evolution or 
the late phase evolution. It should be noted that the majority (70% in a volume- 
limited sample and 77% in a magnitude-limited sample) of SNe la are quite 
homogeneous and only a few objects are true outliers [45| . 

The simple classification scheme of maximum light spectra has to be expanded 
to include various additional parameters (color, light curve shapes, etc.) for a 
complete description of the supernovae. The coming years, with many large 
supernova surveys planned (Sec. will provide the detailed data to define ad- 
ditional classification criteria. 

Several major samples of supernovae have been assembled in the past two 
decades. Leibundgut [46j provided an overview of the situation a decade ago. 
In the meantime, the Center for Astrophysics supernova program has produced 
an impressive sample of around 160 nearby SNe la |47] . the Lick Observatory 
Supernova Search has assembled several hundred SNe la and has published a 
first set recently [IH], the Carnegie Supernova Project has also collected a fairly 
large number of nearby SNe la with both optical and near-infrared data [39]. 
The SNfactory has implemented a bold new technique for the study of super- 
novae: rather than carrying out a series of photometry measurements in indi- 
vidual broadband filters, SNfactory uses an integral field spectrograph to create 
synthesized photometry corresponding to any possible optical filter set. The 
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Palomar Transient Factory [50] and the PanSTARRS supernova searches have 

started and have already produced a number of interesting objects. The distant 

searches wiU be described in Sec. [71 The nearby supernovae are critical for the 

cosmology as they provide the comparison sample and also the anchoring point 

of the distance scale. 

1.2 Supernova physics 

The display of a supernova is the result of the explosion mechanism and the 
structure of the exploding star. Collapsing about one solar mass from the radius 
of the Sun to a neutron star or a black hole releases about 10^^ J, mostly in 
neutrinos. About 10^^ of this energy goes into the acceleration of the stellar 
material and 10~^ into radiation (10^'^ J). The energy gain of burning about one 
solar mass to iron-group elements due to the higher nuclear binding energy also 
is around 10^^ J. Hence the two mechanisms produce very similar luminosities. 

Type II and Type Ib/c SNe are associated with core-collapse supernovae of 
massive stars [51j, while the Type la supernovae most likely are due to the ther- 
monuclear explosion of a compact white dwarf star [52j. 

1.2.1 Thermonuclear explosions The similarity between SNIa has been 
linked to a threshold phenomenon involving at least one white dwarf accreting 
mass from the surrounding environment. These compact remnants of low mass 
stars emit only thermal radiation from their slow gravitational contraction. 

Thermonuclear explosions are due to the explosive burning of carbon and oxy- 
gen. The energy released from the synthesis of iron-group elements is mostly 
used to overcome the binding energy to disintegrate the star. The fairly quick 
brightness evolution (light curve) of Type la supernovae (see Fig. [1]) , a rise time 



10 Goobar Sz Leibundgut 

to maximum light of about 17 to 20 days and a rapid decline after peak indicates 

a small progenitor star. 

The temperatures and densities for carbon burning are reached shortly before 
the Chandrasekhar-mass, the highest mass a non-rotating white dwarf can obtain 
before gravity overcomes the electron degeneracy and collapses the white dwarf. 
Due to the electron degeneracy the mass-radius relation is inverted and the more 
massive white dwarfs actually have smaller radii. Once the carbon burning is 
triggered, it starts a simmering phase within the star, which can last several hun- 
dred years. During this phase the burning creates a large instability, which at 
some point turns into a runaway (e.g. |52]). At this point the burning proceeds 
through the star very rapidly and disrupts it completely within seconds (e.g. [53])- 
There are two ways the burning front can move through the star: at subsonic 
speeds, called a deflagration, or at supersonic velocity, designated detonation. 
Pure deflagrations typically have the problem that they do not achieve enough 
energy for the bright display of Type la SNe, but they have been suggested for 
several of the less luminous SNe la. Detonations have been disfavored for many 
years as they burn most of the star to iron-group elements and do not leave a sig- 
nificant amount of intermediate-mass elements, like silicon, calcium and sulphur 
observed in the spectra. Models with a transition from deflagration to detonation 
have been favored. In this case, the subsonic burning front accelerates to reach 
supersonic speed in the outer layers and hence provides a mix of iron-group and 
intermediate-mass elements. The exact reason for the transition, however, has 
not been completely identified. 

Mass accretion is critical in this scenario. The white dwarf can only grow in 
mass, if it receives matter from a companion star. This immediately requires 
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a close binary stellar system. Also, the white dwarf needs to grow to close to 

the Chandrasekhar mass (~ I.4M0), which is not easy given that typical white 
dwarfs are normally below a solar mass. Possible accretion scenarios, both with 
significant problems at the moment, are either the merging of two white dwarfs 
(e.g. ^54j for at least sub-luminous SNe la) or mass loss from the nearby compan- 
ion star. There are observational arguments for both these channels. A further 
option is to trigger the explosion well below the Chandrasekhar mass by explosive 
burning of surface material. In such a case, helium detonates near the surface, 
possibly where the accretion reaches the white dwarf, and then triggers the central 
carbon explosion by pressure waves focussing near the center. While this possi- 
bility has for a long time been regarded as producing the wrong nucleosynthesis, 
it has recently been revived in full 3-dimensional simulations |55y56j. 

The optical display of a Type la supernova comes from the radioactive de- 
cay of newly synthesized material. Adiabatic cooling of the ejecta means that 
there is no remnant energy for optical emission. Instead the decay chain from 
^^Ni— 7>^^Co— 7>^^Fe (e.g. [57j) is responsible for powering the light curve. The radi- 
ation transport in Type la supernovae is very complicated as the original 7— rays 
from the radioactive decay are down-scattered to lower energies. This happens 
mostly through Compton scattering and then absorption and line emission in 
local atoms (e.g. [58]). This non-thermal process is extremely difficult to model 
and depends on complete atomic data for many ionized atoms. The emission 
at different wavelengths is modulated by these radiation transport effects. A 
striking example is the secondary peak, which is observed in the light curves of 
red and infrared filters (at A > 700nm, see Fig. [T|). This has been tentatively 
explained as due to the ionisation state of iron [59], where the opacities change 
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dramatically for doubly- to singly-ionised iron. 

The original hopes of a single explosion mechanism for all SNe la is probably no 
longer tenable. The observed diversity indicates that several of the evolutionary 
scenarii and also explosion mechanisms are at work. It has also emerged that 
the explosions are unlikely to be spherical and hence viewing angle effects need 
to be taken into account as well |60j . In particular, they can introduce a natural 
scatter into the observables. 

Critical observational input to the models would be the total mass of the ex- 
plosion, the newly synthesized ^^Ni mass and the explosion energy (e.g. [|6l]). 
This information is difficult to come by and is being assembled only slowly. 

1.3 Supernova rates 

Supernovae are rare objects. The last millennium has only six confirmed super- 
novae in our own Milky Way and SN 1987A in the Large Magellanic Cloud. The 
local rate of Type la supernovae is about 3 • 10~^ SNe la Mpc^'^yr^^ [62tl63]. 
The past decade has seen the industrialization of sky monitoring and supernovae 
are now found regularly. The challenge is different for nearby and distant super- 
novae. Due to the limited volume it is rather difficult for find nearby supernovae 
as large sky areas need to be monitored. Nevertheless, the rate of supernova 
discoveries in the local volume (out to about 100 Mpc) has doubled in the past 
decade. For higher redshifts, the increase has been much larger due to the larger 
mirror/longer exposures and wide-area capabilities, since distant supernovae are 
found by deep pencil-beam searches. Due to the great depth, a large volume is 
searched and many supernovae are found. 

The SN la rate is connected to the star formation of the local galaxy and 
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can also provide information on the gestation time before a supernova explosion 

occurs. The progenitor systems of SNe la have not been uniquely identified and 

hence the rate as a function of redshift (or look-back time) provides important 

clues. For systems, which take a long time to grow the white dwarf towards 

the Chandrasekhar mass (several 10^ years) there has to be a cutoff when the 

universe is not old enough to produce a Type la supernova. Figure [2] displays the 

current state of affairs. There is clearly an increase by a factor of 3 in the SN la 

rate back to a redshift of 1. If the drop for z>l is real then we may see the onset 

of the first supernovae at that time. 

Searches at very high redshifts are currently under way (as part of the CLASH 

and CANDELS HST multi-cycle projects). They will show whether indeed there 

are no Type la supernovae at z = 2 (about 3 billion years after the Big Bang for 

the canonical cosmological parameters. 

2 Standard cosmology 

Modern cosmology is based on Einstein's field equations of General Relativity 
(GR): 

R^iu - ^9iiuR = SirCTf^u + ^9(11^, (1) 

which connects spatial curvature (left hand side) and energy content (right hand 
side). Equation [1] includes the extra term A, also called the cosmological constant, 
that Einstein introduced in 1917 to account for the possibility that, on the largest 
scales of the universe, the attractive nature of gravity could be neutralized, thus 
potentially allowing for a static solution. In the context of quantum field theory, 
a A term arises naturally as a vacuum energy density, p^ac = i-e. the 

zero-point energy associated with the production and annihilation of all virtual 
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particles. 



Goobar <fe Leibundgut 



2.1 Priedmann-Lemaitre models 

To compute cosmological observables from GR, Eq. ([1]) is solved using two simpli- 
fying assumptions about the universe: homogeneity and isotropy over very large 
scales. Expressed in the Robertson- Walker line element, distances in space-time 
become (c = 1): 

\1 — kr^ I 
where a{t) is scale factor of the expanding universe and k indicates the curvature: 
(0, 1, —1) for a flat, closed or open universe. 

The solutions to the temporal and spatial components of the Einstein's field 
equations in the Robertson- Walker metric lead to two fundamental equations for 
the dynamics of the universe first found by Friedmann and Lemaitre for the first 
and second derivative of the scale factor, a: 



d\2 SttG n. 

a J 6 



and 



d AttG , 

- = ^(P + 3p), (4) 

a 6 

where p and p are the energy density and pressure of the cosmic fluids, and are 
typically related through the equation of state parameter as: 

p = w • p (5) 

Since the discovery of the expansion of the universe in the late 1920's, the 
current expansion rate is called the Hubble constant, - = Hq. The goal of 



SN cosmology: legacy and future 15 
observational cosmology since then, is to measure the time dependence of the 

expansion rate, H{t), thus indirectly probing the contents of the universe. 



2.2 Components in the Priedmann-Lemaitre cosmological mod- 
els 

In the Standard Model of cosmology, the universe has been expanding and cooling 
since almost 14 billion years (~ Hq^). Thus, after the initial "Big Bang", the 
energy content of the universe was dominated by relativistic particles {w = 1/3), 
eventually surpassed by non-relativistic matter (w = 0) as the temperature fell 
well below the mass of standard model particles^. 

A continuity equation can be derived from Eqs. ([3115]): 

p + 3Hp{l + w) = (6) 

The density of cosmic fluids in the expanding universe evolves according to 
Eq. dlD as: 

p = poe -^"o - , (7) 

where ao is the scale factor today. 

A simple relation is thus found for fluids with constant w: 

P = Po[—) ■ (8) 



a 



We thus find the familiar dilution for non-relativistic matter as the universe 
expands, pM = P% (if) '^5 ^^'^ because of redshift, by another power of a for 



convenient relation between age of the universe and the plasma temperature in the early 
universe is r(MeV) = t~^{s). 



relativistic matter pR = pj^ 
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Since A is constant, is is directly seen from 



Eq. ([8]) that its equation of state is wa = — 1. More generically, any dominant fluid 
with w < —1/3 will cause the universe to accelerate, as seen from Eq. (jl]). Such 
fluids are called dark energy, the "simplest" one being the vacuum energy density 
associated with the cosmological constant, A. Other dark energy candidates are 
discussed in e.g. |5t[T3]. 

It is customary to re-write the Friedmann equation ^ in terms of the cos- 
mological redshift, 1 + z = and the density terms normalized by the critical 
density corresponding to a flat universe, ^Ide = ^^jf- = Pde ■ and the cur- 
vature term is defined as Qk = —k{aoHo)~'^ . 

Since we lack the knowledge of the nature and properties of the equation of 
state of dark energy, phenomenological parameterizations of w{z) are used for 
comparisons with data, e.g. 



where the free parameters wq and Wa are fitted with observational data. 

2.3 The luminosity distance and the Hubble diagram 

There are several cosmological distances measures. Here we concentrate on the 
Luminosity Distance, di, which is defined from the expected flux, J^, from an 
object with luminosity L (in units of energy per unity time and unit area), = 
Since the universe expands, any cosmological distance estimate will depend 
on the growth of the scale factor and thus on the Hubble parameter. For a full 
derivation, the reader should consult a cosmology text book, e.g. |68) . Here 
we concentrate on the final expression (neglecting the radiation density term 




(9) 



10-5): 
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(1 + z)-7^sin(7^rT^I), QkKO 

dL = { {l + z)I, nK = o (10) 

(1 + z) 1 sinh(Vn]F/), nK>o 

Qk = i-nM-^DE, (11) 



Hiz) = Ho.J{i + z)^nM + fiz)nnE + {i + zynK, (i3) 

,i + w{z'y 



f{z) = exp 



3 / dz 



(14) 



/o 1 + ^' 

It is customary to use logarithmic distances in astronomy, the distance modulus, 
H = 51og;^o '^i(Mpc) + 25 = nif — Mf, where the last term corresponds to the 
difference between the measured magnitude of the object, ruf and its absolute 
magnitude through the same filter /, Mf, i.e. corresponding to the flux of the 
object if placed 10 pc away from the observer. Cosmological model testing is done 
by comparing measured distance moduli vs. redshift against model predictions. 
This is usually called the Hubble diagram. The supernova Hubble diagram, which 
will be discussed in more detail in the following is displayed in Fig. [31 

3 Cosmology with Type la supernovae 

Since SNe la display rather uniform light curves, they can be used as local clocks. 
This test for cosmological time dilation had been proposed soon after supernovae 
were recognized and was implemented with the first distant supernovae |701l71j. 
The spectral evolution also tests for this stretching of time in the observer frame 
|72j . All tests so far have confirmed the prediction of general relativity. The 
corrections for time dilation are now routinely applied to the supernova data. 
Most supernova searches operate by repeated imaging of patches of the sky 
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with many thousand galaxies. Supernovae are found through pixel-by-pixel sub- 
tractions of CCD images obtained at different epochs, as shown in Figure HI 

Once a supernova is detected, a few steps are necessary before the brightness 
can be used for distances estimates. First and foremost, the object is spectro- 
scopically classified, as described in Sec. II. 11 The spectrum also provides the vital 
redshift for the object. In many cases light from the host galaxy can be analyzed 
to find a redshift, but for "hostless" supernovae (where the surface brightness of 
the host galaxy falls below the detection limit) the redshift can also be determined 
from the supernova spectrum itself. Figure [5] shows the identification spectra of 
two Type la supernovae, SN2003du in the nearby universe, and SN2000fr at 
z = 0.543, i.e., more than 8 Gyrs apart. The striking spectral similarity is the 
key of the "standard candle" technique. Since distances are established based 
on the peak flux, light curves are built through several filters, ideally including 
observations well before and after maximum, as shown in Fig. [T] Besides time 
dilation making the light curve width increasingly broader with redshift, SNIa 
optical light curves exhibit a range in widths. The correlation between light curve 
and supernova luminosity has been established over several years. More luminous 
supernovae display typically a slightly slower evolution. 

Intervening dust can scatter and absorb supernova light and hence produce a 
dimming, which needs to be corrected for the distance measurement. Additional 
effects may also infiuence the observed colors of SNIa, as discussed in Section [5l 
Several light curve "fitters" have been proposed to extract the key parameters 
from the data: the day of maximum and the corresponding peak flux, the light 
curve shape, and the reddening (color excess). The latter is often computed as the 
relative magnitude attenuation between the restframe i?-band and the V-band. 
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The empirical color-brightness and lighcturve shape-brightness relations used to 

"standardize" SNe la are shown in Fig. [6l 

Since the measurements are done with fixed filter pass-bands, irrespectively 
of the supernova redshift, the measured fluxes correspond to diff'erent parts of 
the restframe spectrum. The distance estimates rely on flux ratios of the same 
restframe fllters and hence, a ET-correction is needed to transfer the observed flux 
into the restframe for comparison with local objects. Due to the signiflcant spec- 
tral evolution, high S/N spectroscopic observations would be needed for every 
flux point, a daunting task! However, given the uniformity of SNe la, spectral 
templates are used to compute X-corrections. These templates are made of well 
sampled, low noise, measurements of SNe la, mostly at low redshift. Any poten- 
tial reddening is determined by comparison of the multi-band data with synthetic 
colors calculated from spectral templates after applying an extinction law, e.g., 
the one found for stars in the Milky Way [7H]. The photometric measurements 
from different telescopes also need to be corrected for differences in the wave- 
length response of the telescope reflectivity and instrument sensitivity. Modern 
analysis methods try to correct for these systematic effects in the light curve 
fitting. 

To summarize, observed SN peak magnitudes at redshift z through a filter Y 
are compared with model predictions taking into account the observed color and 
light curve shape (c, s). The effect of cosmological redshift on the spectrum, from 
the restframe wavelength X to the observations in Y are handled through the 
cross-filter Xxy-corrections [77]. These would typically also include effects from 
reddening and light curve shape. 
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niYiz, d, c) = Mx{s) + 5log^o dUz) + KxYiz,s,c) + A^^ , (15) 

where the last term corresponds to dimming by Milky- Way dust. Note that in 
Eq. (jlSp it is implicitly assumed that corrections for extinction in the host galaxy 
are included in Kxy- 

3.1 The Hubble constant 



As discussed in Sec. 12.31 the relation between the distance modulus and cosmo- 
logical redshift can be used to measure cosmological parameters such as Hq, 0,m, 
i^DE and the dark energy equation of state parameter, w. For z <C 1, we find 
from Eq. (fTi]) : 

dL = ^^(l + l{l-qo)z + 0{z')y (16) 

where we have explicitly reinserted the speed of light c into the expression and 
introduced the deceleration parameter, qo, which is related to the cosmological 
"fluids" as: 

The linear part of the Hubble diagram, cIl^z <c 1) ~ cz/Hq, has been used to 
measure Hq since several decades using magnitude measurements to SNIa using 
optical data, e.g., through the i?-band filter, uib' 

^ = me - Mb = 5 log cz - 5 log Ho + 25 + Kbb + ^. (19) 

Thus, log Hq can be read off from the intercept of the log cz — 0.2mB relation (cf. 
Fig. El). 
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While SNe la are very useful for measuring relative distances, extracting Hq 

requires an absolute calibration of the SNIa brightness. Since theoretical models 

are not sufficiently accurate, the SNIa absolute magnitude needs to be calibrated 

observationally [22j . Cepheid stars are used to measure distances to galaxies in 

the nearby universe. By now, Cepheid distances to 15 SNIa host galaxies within 

about 20 Mpc have been obtained with the Hubble Space Telescope [78]. Riess 

et al. [79] mitigated some of the shortcomings in building up the distance ladder 

to obtain absolute distances to SN la hosts. Near-IR observations of Cepheids 

in six systems with accurate SN la data were used to improve the calibration 

and reducing the uncertainty from extinction by dust. Furthermore, a better 

anchoring system was found to calibrate the Cepheid distances, the "maser" 

galaxy NGC 4258, for which an accurate absolute distance can be derived. Riess 

et al. [79] report Hq = 74.2 it 3.6 km s~"'^Mpc~^, where the 4.8 % error includes 

systematic uncertainties. 

3.2 The accelerated expansion 

The extension of the SNIa Hubble diagram to higher redshifts, needed to probe 
the change in the cosmic expansion rate, requires a specific observational strategy: 
large patches of the sky need to be scanned regularly. This prompted the need for 
medium size telescopes, 2 — 4 meters in diameter, equipped with CCD cameras 
with large fields of view. 

By the mid-nineties, Goobar & Perlmutter [80] showed that constraining the 
value of the cosmological constant. A, was feasible using SNIa data, if the surveys 
targeted a wide range of redshifts, thereby allowing to break the degeneracy in 
cosmological parameters. 
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Figure [7] shows the il,M — ^A plane, with the two original results of the SCP [16] 

and HZ SN Team [15j, along with the recent Union2 compilation based on 557 

SNIa I69j. The data-set includes 192 SNe with z < 0.15 from gTlEIHHl]. The 

intermediate and high-z includes samples from |15M16|l73ll85H92j . adding up to 106 

objects with 0.15 < z < 0.55, 64 with 0.55 < z < 0.75, 41 with 0.75 < z < 0.95 

and 28 SNIa with z > 0.95. The SN contours are derived from the Hubble 

diagram data shown in Figure [3l 

4 Cross-cutting techniques 

The power of the SNIa data for constraining cosmological parameters is vastly en- 
hanced when the likelihood functions are combined with what is found with other 
probes, most notably the cosmic microwave background (CMB) anisotropics and 
baryon acoustic oscillations (BAO), as shown in Figure [7] and discussed in [12] . 
The combination of the different techniques becomes dramatically clear in the 
determination of w. For this measurement all methods show large degeneracies, 
but their combination provides a powerful tool to constrain a constant w, as 
shown in Fig. [8l 

Similarly to the use of SNIa as "standard candles" for distance estimates, the 
CMB and BAO techniques rely on "standard rulers". Matter-radiation oscilla- 
tions with wavelengths given by the product of the sound speed and the elapsed 
time since the Big Bang leave a characteristic bump in the power spectrum of 
the CMB temperature fluctuations originating at z ~ 1100 and corresponding to 
an angular scale of about 1° on the sky today. Measurement of this angular scale 
in WMAP 7- year data [9lj leads to the precise constraints cosmological parame- 
ters, as shown in Figs. [7] and [51 As the universe expands and photons decouple 
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from baryons, the matter anisotropies grow into the structures observed today, 

down to the scale of galaxies and clusters of galaxies. By studying a sample of 

Luminous Red Galaxies in the SDSS survey, Eisenstein et al. |95j could measure 

the ~ 150 Mpc length of the "ruler" at 2: ~ 0.35 leading to significant parameter 

constraints, mainly on VLm- 

Whereas the SNIa Hubble diagram is anchored at z ~ 0, the CMB and BAO 

originate in the early universe. As a result, the "orientation" of the confidence 

regions differ between the techniques, thus rendering the crossing of the allowed 

regions particularly useful to derive accurate cosmological parameters. 

5 The era of precision cosmology 

The pioneering efforts by the SCP (led by Saul Perlmutter) and the High-Z team 
(led by Brian Schmidt) in the mid and late 90's were followed by new SN programs 
which developed the field in several directions. 

Two major programs were established to collect large samples of high quality 
data in the z = [0.3— 1] redshift range: the ESSENCE program at the CTIO 4-m 
in Chile [96] and SNLS at the 3.6-m CFHT in Hawaii [89]. These programs were 
backed up with significant spectroscopic identification efforts from 8 — 10 meter 
class telescopes: VLT, Keck, Gemini and Magellan. The 3-year SNLS sample 
recently reported by [97j includes over 250 SN la, i.e., about 180 additional su- 
pernovae to what is already included in the Union2 sample. SNLS and ESSENCE 
have additional ~300 SNe la unpublished so far. 

One of the observational breakthroughs of SNLS was introducing the "rolling 
search" technique. One-square-degree fields were revisited every 4 to 5 nights 
in multiple bands, thus new objects were found while light curves were built of 
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previously discovered supernovae. The excellent time sampling with several filters 

and with a single telescope and CCD system dramatically improved the quality 

of the SN light curves. 

The redshift "gap" around z ~ 0.2 is filled by the SDSS-II supernova survey. 

The first data [98j includes 103 SNe la with redshifts 0.04 < z < 0.42. Overah, 

it is expected that the SDSS-II supernova survey has about 400 supernovae for 

the analysis. 

Finally, the Hubble diagram is anchored with supernovae from the Calan-Tololo 
[HT] and CfA surveys [171[82l[83] (see also the list of nearby projects given in 
Sec. [LID. 

The HST program led by Adam Riess j90tl99j dramatically enlarged the ob- 
served redshift range of SNe la. Thanks to the near-IR capabilities from space, 
the z > 1 domain was opened, allowing the study of SNe in the decelerated phase 
of the expansion, i.e. exploring the transition between acceleration and decelera- 
tion marked by the turn-around of the differential Hubble diagram in the bottom 
panel of Figure [3l 

Progress at z > 0.9 is slower, mainly because it hinges on major access to 
near-IR observations with HST. The SN detection efficiency can be improved by 
targeting massive galaxy clusters, as demonstrated by |100j . Currently, a massive 
project to find very distant {z > 1.5) supernovae has started making use of two 
multi-cycle HST projects, which image several areas and galaxy clusters, over 
several years. Yet, it is unlikely that the number of z > 1 SNe la that can be 
used in the Hubble diagram will exceed ~ 100 with current instruments. 

The cosmological results are summarized in Table [TJ 

The results are compatible and the strongest deviations come from the use 
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of the light curve fitters (and their assumptions on reddening; see Sec. 16. 4p . It 

should also be pointed out that the samples are not completely independent. The 

nearby supernova sample for almost all results presented in Tab. [T]are identical. 

The Union2 data set is the collection of all published SN la data suitable for 

cosmology. 

6 Systematic uncertainties 

6.1 Calibration 

Precise distance estimates with Type la supernovae require very accurate relative 
photometry (Aw ~ 2 Am). While future surveys are likely to measure most (all) 
their SNe with a single, well calibrated system, current cosmological bounds are 
derived from a compilation of surveys. There is a worry for calibration offsets 
between the data sets, potentially resulting in a bias in the fitted cosmology. The 
SNLS team has been extremely active in refining the calibration of the CFHT 
MegaCam system and have reached sub-percent precision in their photometric 
calibration [101]. Unfortunately, significant calibration uncertainties remain for 
the near-IR HST system, including NICMOS nonlinearities. This amounts to 
a sizable contribution to the budget of systematic errors for the z^l SNe la 
measured up to now |102| . Clearly, future SN la programs must integrate a 
robust calibration program into the survey design. 

6.2 The UV-flux 

Since high-z surveys are mostly confined to optical bands, the measured (red- 
shifted) fluxes originate at ultraviolet and blue wavelengths. For restframe wave- 
lengths shorter than what corresponds to the S-band, several sources of concern 
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have been raised. The atmospheric cut-off and steeply falling instrument sensitiv- 
ity has lead to significant difficulties in providing a good estimate of the intrinsic 
brightness in the U -band. About 5% of the SNIa flux is emitted at wavelengths 
shorter than the [/-band, not at all accessible from the ground for z = 0. The 
wavelength region below 3000 A remains one of the least explored, in spite of the 
fact that it is rich in key information, e.g., about the underlying explosion physics. 
Furthermore, since the dimming by dust grains increases rapidly at shorter wave- 
lengths, the uncertainty in host galaxy dust extinction properties and intrinsic 
color variations are major limitations for precise distance measurements. 

Photometric studies with the Swift satellite |1031ll04j and STIS spectroscopy 
with HST |105| have been used to test model predictions suggesting an increased 
dispersion in the UV, because of its larger sensitivity to progenitor composition 
|106H108| . These studies, as well as similar ones at z ~ 0.5 |109| do indeed find 
a scatter about twice as large in the UV compared to the optical. 

It remains unclear what these effects mean for the use of SNIa as distance 
indicators. 

6.3 Reddening and absorption 

Accurate corrections for absorption of the SN light along the path from the source 
to the telescope are crucial for distance determinations. Interactions of light and 
dust grains may happen in the circumstellar material around the exploding star, 
in the interstellar medium of the host galaxy, in the intergalactic medium and 
in the Milky Way. Absorption in the Earth's atmosphere is corrected through 
the photometric calibration. The distribution and properties of dust grains in 
the Milky Way is known well enough to allow for correction accurate at the 
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sub-percent level. The other terms remain more problematic. The dust size- 
distribution is captured by the Galactic extinction law and is often described by 
the total-to-selective extinction ratio, Ry- For the Milky Way, an average value 
Rv = 3.1 has been found |76j . with relatively small scatter around this value 
across the galaxy. The existence of significant amounts of intergalactic dust, 
^dust ~ 10~^ — 10~^, has been inferred from estimates of the stellar density and 
metallicity as a function of redshift. After the evidence for cosmic acceleration 
was presented in 1998, several authors studied the possibility that the faintness of 
SNe la at high-z could instead be caused by "gray" interstellar dust in the inter- 
galactic medium |110H112] . An upper limit on the restframe i?-band attenuation 
^^(z = 1) < 0.1 on impact of dimming by dust in the intergalactic medium was 
derived by studies of quasar colors |113yil4j . 

More recently, Menard et al. |115j performed a cross-correlation study between 
colors of nearby and distant quasars and reported a statistical detection of dust 
reddening up to large scales around galaxies. From this result, they inferred the 
opacity of the universe. Their model-dependent estimate gives Ab{z = 0.3) ~ 
0.02 mag. This attenuation is partially compensated by the color-brightness 
correction. However, since the attenuation happens along the line of sight rather 
than at the redshift of the SN, the current light curve fitters cannot account for 
this effect in its totality as they correct for discrete redshifts of the host galaxy. 

Considerably larger extinction by dust may occur in the interstellar medium 
of the host galaxy, especially for SN explosions close to the galaxy center. Multi- 
band photometry of SNe la has been used to constrain the wavelength dependence 
of the absorption. There is by now a large number of studies reporting measure- 
ments of Ry ~ 2 based on SN la colors, (e.g. |116H121] ). Moreover, low- values 
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of Rv do also provide a smaller scatter in the Hubble diagram and thus a better 

fit in global parameter fitting of SNe la across a wide redshift range. This is 

somewhat puzzling and indicates that reddening by dust is not fully understood. 

Multiple scattering of photons has been suggested as possible explanation to 

the "unexpected" reddening law. This could be the case in the presence of dust 

in the circumstellar environment of the SN la progenitor, as suggested by [122| . 

who also provided an effective reddening law that gives good fits to the best 

measured/most reddened SNe la in the CSP low-z sample [l9]. Spectroscopic 

evidence for circumstellar material around highly reddened SNe la has been found 

in e.g. [T23] . 

At this point, the lack of understanding of the relative strengths and specific 
properties of the various effects contributing to the color-brightness relation of 
SNe la constitutes one of the biggest concerns, especially since a redshift depen- 
dence cannot be excluded. 

6.4 Light curve fitters 

Establishing the peak brightness of SNe la has to date been essential essential 
for their use in cosmology. A major progress over the past decade has been 
the much more systematic sampling of supernova observations. While the light 
curves of early supernovae had to be established from rather sparsely sampled 
data, the current and future surveys provide a much denser temporal coverage 
and the determination of the peak brightness is highly improved. However, light 
curve shape and reddening have to be established through fitting. Since color and 
shape are measured with respect to some reference set and may also include some 
additional priors, several light curve fitters have been proposed over the years. 
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Examples are the fitter developed by the SCP, SNMinuit, used the original light 

curve shapes in [124J based on low-z data, and improved in [71j. A "stretch" 

factor, s, was introduced to characterize the light curve shape perturbations with 

a simple scaling of the time axis within the, t ■ s, and a corresponding brightness 

correction AM = a ■ {1 — s), i.e. broader light curves s > 1 corresponding to 

brighter intrinsic brightness, AM < (cf. Fig. [6|). While s is fitted with points 

both in the rising and falling parts of the light curves, other shape estimates, 

like Amis |125j . use only a limited region after maximum. The two major LC 

fitters used today are SALT2 |118) and MLCS2k2 |126j . the second generations of 

SALT [127j and MLCS |128j . The major differences between the two reside in the 

reference spectral templates in the blue and UV regions and in the interpretation 

of reddening. While the MLCS approach is to obtain the reference samples 

from the low-z SNe, SALT combines data from both high and low redshift. In 

consequence, the MLCS "training" may be biased by the difficulties to obtain 

high quality data in the UV, as explained in section [6^ whereas SALT may lack 

the ability to detect intrinsic changes with redshift and makes use of faint data. 

Another key difference involves the treatment of color-brightness relation. MLCS 

assumes that SN la reddening is caused by dust, thus negative color corrections 

are not used. In contrast, SALT combines all possible contributions into a single 

function relating the intensity as a function of wavelength, and both positive 

as well as negative corrections are allowed, i.e. a negative color excess implies 

increased brightness and reddening corresponds to dimming of light. An example 

of MLCS2k2 and SALT2 fits of the SNLS supernova 04D2gp {z = 0.732) is 

shown in Figure [9l Although both fits give good x^i SALT2 finds the SNe to be 

moderately blue, while MLCS2k2 puts it at small (but positive) extinction. Thus, 
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the derived distance modulus for this particular SN differs by 0.1 mag [98]. More 

recently, two other fitters have been developed SiFTO [129J, so far only used by 

the SNLS collaboration, and SNooPy [64j by CSP that, unlike other fitters, can 

also be used to build near-IR light curves. 

6.5 Selection bias 

Varying selection efficiency is an inherent difficulty in a technique involving com- 
parisons of brightness of objects at different redshifts. In particular, cosmological 
estimates using multiple surveys need to model (most often through Monte Carlo 
simulations) the fraction of the SN la luminosity function sampled by the de- 
tection requirements of each instrument. At high-z, the "bottle-neck" is the 
spectroscopic identification, which is only practically obtainable for i < 25 mag 
in 8-10 meter class telescopes. This naturally leads to a redshift dependent selec- 
tion bias in magnitude and color, e.g., for the 3-year SNLS sample [97], as shown 
in Figure [TOl Perrett et al. |130j examine these selection effects in detail, and the 
related corrections. A problem associated with these corrections is our lack of 
precise understanding of potential drifts in either reddening properties (Sec. I6.3p 
or intrinsic SN la luminosity, to be discussed further in Sec. 16. 8[ An ambitious 
space mission equipped with a multi-band imager and a sensitive spectrograph 
could potentially resolve this difficulty (see Sec. [8|). 

6.6 Gravitational lensing 

In Section 12.11 we derived the expression for the luminosity distance used to 
measure cosmological parameters assuming that the universe is homogeneous 
and isotropic. Although there is ample empirical evidence that this is a good 
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approximation at the largest scales, e.g., from the cosmic microwave background, 

it is clear that inhomogeneities exist at smaller scales. Thus, foreground galaxy 

halos and clusters of galaxies act as gravitational lenses and introduce a "scatter" 

in the observed SN la flux that increases with redshift |13H - I133] . Since the 

number of photons is conserved under lensing (some objects get magnified and 

others de-magnified), the luminosity distance averaged over a large number of 

SNe la is expected to converge towards Eq. (jlOp . However, for small number 

statistics or if the sample is magnitude limited, gravitational lensing could be a 

source of selection bias. This was first investigated in |134j for the HST GOODS 

sample ^99j, where the redshifts of the foreground galaxies are well measured 

and their masses can be estimated through their luminosities and the effect of 

lensing of individual SN can be evaluated. More recently, Jonsson et al. [135] have 

studied the weak lensing perturbations in the SNLS 3-year data (Figure [TT|) . The 

distribution of lensing magnifications of SNe la as a function of redshift matches 

the expectations and no significant bias has been detected in the current high-z 

sample. 

6.7 Peculiar velocities 

Local matter inhomogeneities are also important for precision cosmology with 
SNe la. Local density perturbations induce peculiar velocities in galaxies of 
about 300 km s^^. In addition, the exploding star could have a comparable 
internal motion within the host galaxy. Since these velocity components cannot 
be disentangled from the cosmological redshift, an extra source of uncertainty 
needs to be considered, cr^ ^ • jj^^^ 
intrinsic brightness cr^ ~ 0.1 mag for z < 0.015. The concern, however, is not 



l+z 
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which exceeds the error from 
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the addition of an added "random" error, but rather the possibility of an effect 

that does not average out with statistics. In particular, our own motion and 

the correlated peculiar velocities of the SN la host galaxies at low-z could have 

a non- negligible effect on the derived cosmological parameters. Several authors 

have studied these effects, concluding that inclusion of SNe la down to z = 0.02 

may result in a systematic uncertainty of about Aw^O.Ol [136J. 

6.8 Brightness evolution 

The brightness distribution of SNe la has been shown to depend on the host 
galaxy environment. E.g. Hamuy et al. [137j noticed that the average B and 
V-band light curve width differs between galaxy types. It is now an established 
fact that elliptical galaxies host a larger fraction of narrower, low "stretch", su- 
pernovae. Furthermore, recent studies of nearby supernovae [138) : SNLS super- 
novae |139 j and SDSS SNe |140j suggest possible evidence for varying supernova 
properties with host galaxy not entirely corrected by the standard width-color- 
brightness relations. In particular, it has been argued that splitting the sample 
according to host-galaxy mass and allowing the peak absolute magnitudes of the 
two samples to differ improves the residuals in a statistically significant manner. 
The evidence suggests that large stellar mass galaxies host the brightest SNe la 
(AM ~ 0.075 mag), after light curve shape corrections. Given the time evolu- 
tion of metallicity and stellar mass of galaxies, the possibility of a drift in the 
brightness of the "standard candle" with redshift remains a source of concern for 
precision cosmology. 
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6.9 The "averaging" uncertainty 

After the discovery of dark energy several authors |141H145] raised doubts about 
the theoretical foundations leading to Eq. (jl4|) . In particular, it has been claimed 
that the average expansion in a locally inhomogeneous universe behaves differ- 
ently than expected in Friedmann-Lemaitre models based on the large scale aver- 
age energy density. The reason being the non-linearity of the Einstein equations: 
spatial averaging and solving the Einstein equations do not commute |146] . This 
is sometimes referred to as the "back-reaction" problem. Although the issue 
is generally accepted, arguments have been put forward for why the effect is 
too small to explain "cosmic acceleration" without dark energy. In a matter- 
dominated universe the local dynamics should be (almost) Newtonian, as long as 
potential perturbations and peculiar velocities are non-relativistic. Since Newto- 
nian gravity is linear, averaging and evolution do commute in the Newtonian limit 
and should commute to a good approximation in general relativity |147M148 j . Any 
relativistic corrections should be much to small to mimic effective acceleration. 
Since it is notoriously difficult to estimate the size of the effects from non-linear 
back-reaction in general relativity, it remains unclear at this point whether this 
presents a limitation for precision tests of w. 

7 Recent developments 

7.1 Expanding the wavelength window 

The use of Type la supernovae helped establish a major milestone in precision 
cosmology, especially the discovery of dark energy. However, further progress in 
the field is limited by our ability to resolve the systematic uncertainties outlined in 
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Section [6l Some outstanding issues involve understanding the progenitor system 

leading to Type la supernovae and the properties of dust extinction affecting 
SN la. A natural path that has been followed for this purpose is to expand the 
studied wavelength window, especially for nearby SNe. Furthermore, theoretical 
arguments have been put forward suggesting that the dispersion of intrinsic SNIa 
peak magnitudes should have a minimum around 1.6;um |59) . 
7.1.1 The near-IR The Earth's atmosphere severely limits the feasibility 
of ground based observations in the near-IR. Vibrational bands of H2O and CO2 
block the incoming radiation beyond 1 ;um, leaving a few observational windows 
in the near-IR: J,H and X-bands, centered at 1.2, 1.6 and 2.2 /xm, respectively. 
To make matters worse, significant thermal atmospheric emission (mostly from 
multiple narrow OH-lines) reduces the obtainable signal-to-noise ratio. This is 
particularly challenging since the fractional flux in the near-IR is small for SNIa: 
2, 1 and 0.5 % for J, H, K [H]. In spite of these difficulties, significant efforts have 
been devoted recently to studies of supernovae in the near-IR. A non-standard 
extinction law has been measured with great accuracy for a number of highly 
reddened SNIa with high S/N near-IR data [MKIIZKIISKmKIlg]. Moreover, 
statistical tests on z < 0.05 SNIa in [39 tll50l - ri52j offer observational support for 
the claim in [59j of a somewhat narrower intrinsic scatter in the near-IR. However, 
this benefit may only be exploited at high-redshift with space based observations, 
e.g. the James Webb Space Telescope. Because of the paucity of data in near- 
IR, and the very narrow redshift range studied, it is yet unclear if problems 
with i^-corrections may arise, i.e., varying spectral features moving in and out 
of the filter pass-band as the redshift increases. An interesting development 
taking place in recent years, and more easily accessible than the restframe NIR, 
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is building a Hubble diagram in restframe /-band, i.e., at about 8000 A |153pi54] . 

In particular, Freedman et al. |155j use a sample of 56 SNIa; 21 nearby, and 35 in 

the redshift range 0.1 < z < 0.7, and find limits on w which are consistent with 

the results at shorter wavelengths. 

7.2 Spectral indicators as secondary calibrators 

Thanks to the availability of large spectroscopic sets of SNIa, both nearby and at 
cosmological distances, data mining efforts have been conducted trying to identify 
spectral indicators capable of refining the intrinsic dispersion of the "standard 
candle", e.g., [33 |IMl[T09tll56H168j . Although some of these findings could have 
a profound impact in our understanding of the physics behind SNIa explosions 
and optimize their use in cosmology, the statistical evidence and interpretation 
of the correlations needs further assessment. 

8 Supernova cosmology for the next decade 

The goal of the next supernova cosmology surveys must be to reach a level of 
accuracy promising to meaningfully constrain a time-dependent w{t). From the 
descriptions of the systematics of the supernova cosmology, it is clear that a strict 
control of the surveys is required. This has led to the concentration on singular 
telescopes and their calibration for the searches and the photometry. Finding 
potential further correlations in the SN la observations can further sharpen the 
acuity of SNe la for cosmological distances and will also give hints on the super- 
nova physics. Another improvement would be to explore the small luminosity 
scatter of SNe la in the near-infrared. Building a well-sampled Hubble diagram 
at near-infrared wavelengths would reduced the uncertainties on the reddening 
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corrections as well. 

Several programs aiming at discovering thousands of SNe la up to 2; = 1 are in 
the process of starting up. The "rolling search" technique, where CCD cameras 
with very large field of view are used to revisit sky pointings with cadences of 
a few days, has been proved very successful, e.g., by SNLS using the 1-square- 
degree at the 3.6m CFHT in Hawaii. New SNe are found while multi-band light 
curves are built for transients discovered in earlier epochs. With the imminent 
start of operations of the Dark Energy Survey (DES) on the CTIO 4m (Chile) 
using a new 3 square degree mosaic CCD camera, about 3000 SNe la are expected 
within 5 years. The redshift range will overlap with the existing samples from 
SDSS, SNLS and ESSENCE. 

A fully developed Pan-STARRS with four 1.8m telescopes observing simulta- 
neously will have the capability to scan the entire sky within just a week. With 
the 3-square-degree CCD array on each unit and a pixel size of 0.3 arc-seconds, 
a large number of SNe la can be found across a wide redshift range. The first 
prototype unit, Pan-STARRS 1, is already collecting data and several supernova 
discoveries have been reported. 

The 8-meter telescope foreseen in the LSST project, carrying a 9-square-degree 
camera and performing repeat all sky imaging in five optical bands will propel 
the art of finding supernovae to a whole new dimension. With an expected yearly 
yield of 250000 SNe la, with light curves sampled every five days, low- and high- 2; 
supernovae will be found with the same instrument and calibrated to about 1% 
photometric precision. 

The limitation of this project lies also in the big numbers: it is impossible to 
get spectroscopic identification for more than a small subset of supernovae. The 
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redshift determination most likely will be done through photometric redshifts of 

the host galaxies. Thus, the success of the LSST SN la program for cosmology 
depends on the accuracy of the photometric identification and control of system- 
atic uncertainties related to e.g., a possible evolution of the standard candle and 
brightness-color corrections, as discussed in Secl6j 

First proposed in 1999, the SuperNova Acceleration Probe (SNAP) was for a 
long time the leading concept for a dedicated satellite mission aiming at the study 
of the properties of dark energy with SNIa. Several competing designs have been 
proposed in the mean time, focusing on one or several cosmological probes: Type 
la supernovae, weak lensing measurements and baryon acoustic oscillations. At 
present, two missions are being considered: the European- led EUCLID satellite 
and a US-proposed WFIRST. ESA's EUCLID is mainly focusing on baryonic 
acoustic oscillations and weak lensing and the current specifications include only 
a broad optical ("visible") and three near- infrared filters, which is arguably non- 
optimal for precision studies of SN la over a wide range of redshifts. The WFIRST 
mission concept includes SN la (as well as baryonic acoustic oscillations and 
weak lensing) and a 1.5- meter class telescope has been proposed as a reference. 
Although a space mission cannot compete with LSST in statistics of SNe la, 
it is likely to collect a much better studied sample of SN la, including near-IR 
observations and possibly spectrophotometry, thereby having a better control of 
systematic uncertainties, the ultimate limitation in studying the properties of 
dark energy with SNIa. 

Figure [12] shows the projected sensitivity of WFIRST in the (wo^Wa) param- 
eter space compared with the current constraints from the Union2 sample in 
Amanullah et al. (2010) [691. The shaded "petals" in the figure show classes 
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of "quintessence" models, i.e., where a scalar field dominates the current en- 
ergy density of the universe One such case is the pseudo Nambu-Goldstone 
boson (PNGB) |169j . e.g., an axion with a Hubble scale (Hq^) Compton wave- 
length, i.e., a mass m ~ 10"^'^ eV. Predictions from String Theory inspired scalar 
field potentials in gravitational theories with extra dimensions (SUGRA) and the 
braneworld model by Dvali and Turner [170j are also shown for reference. Clearly, 
any deviation from A would be dramatic breakthrough in fundamental physics. 
Observational cosmology, and distance estimates with SNe la in particular, will 
continue to contribute to the full picture of high-energy and particle physics. 
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CSP SN2007af 




20 40 60 
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Figure 1: Optical and near- infrared light curves of the nearby SN 2007af obtained 
by the Carnegie Supernova Project (Stritzinger et al. in prep.). The light curves 
trace the brightness evolution in six different wavelength passbands ranging from 
~3500A (u) to 1.6 (H). For clarity the individual filter light curves have been 
offset by the amounts indicated. Light curve fits }64j are shown and the key fit 
parameters (time and magnitude at maximum, the width of the light curve and 
a typical color "B-V") are marked. The development of the second peak with 
increasing wavelengths is very prominent. Note that the light curve widths are 
used for comparisons within individual filters. 
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Figure 2: The evolution of the Type la supernova rate per unit volume as 
a function of redshift (data from [65j). The error bars reflect the statistical 
uncertainties per redshift bin. The highest redshift bin contains 3 objects. The 
increased rate with redshift follows the increased star formation rate observed in 
the universe. The two models |661l67j show the effect of different progenitor life 
times. A long progenitor evolution towards a supernova (dash-dotted line) means 
that no SNe la are observed in the early universe. In the case of a shorter delay 
time (dotted line) the rate remains high at high redshifts (and large look-back 
times). 
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Union2 data-set 




Redshift 



Figure 3: The Hubble diagram of Type la supernovae correlating distance 
modulus (/i) vs. redshift. The Union2 compilation |69j represents the currently 
largest SN la sample. The linear expansion in the local universe can be traced 
out to z<0.1. The Hubble constant sets the absolute level of the data, while 
it is irrelevant for the determination of JIa, which is a relative measurement 
(bottom panel). The distance relative to an empty universe model {fiempty] 
Qm = = 0) is shown in the lower panel. The data are binned for clarity 
in this diagram. The blue curve shows the expectation from the best fit ACDM 
model {Qm = 0.3, JIa = 0.7). 
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Figure 4: Detection principle for supernovae: images separated by between a few 
days and weeks are subtracted from each other (NEW - REF = SUB). In this 
case, 5(-band images used to detect SN 1999du {z = 0.260) in [73] are shown. 
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Kessler et al. (2009) 
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Amanullah et al. (2010) 
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^-0.016 


n QQ7+0-050 +0.077 
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Table 1: Cosmological parameters from the first year papers of SNLS, ESSENCE 
and SDSS and the Union2 sample, Refs. [691 [891 [961 [98]. The fits include lower 
redshift data as well as BAO and CMB information. A flat w—CDM model with 
constant w is assumed in all analysis (see text for more detail). 
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Figure 5: The observed spectrum of SN 2003du [71] compared with the restframe 
spectrum of SN 2000fr (z=0.543; |75j). both about 6 restframe days prior to B- 
filter maximum. The two spectra display remarkable similarities in the absorption 
and emission regions. Currently available data have not provided any evidence 
of clearly different spectral appearance of distant SNe la compared to their local 
counterparts indicating that there are no major changes in the explosion physics. 
The spectrum of the nearby supernova extends to much longer wavelengths, which 
are redshifted outside the observable window for the distant supernova, where the 
restframe UV region (not fully shown) is covered much more extensively. The 
bottom panel indicates the optical filter bands for the Johnson-Cousins UBVR 
system as well as the SDSS ugri filter set. 
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Figure 6: The color vs. luminosity (upper panel) and the light curve shape vs. 
luminosity (lower panel, after color corrections) correlations of Type la super- 
novae. Light curve fits to the restframe B filter observations of 685 SNe la |69j 
covering a redshift range from 0.025<z<1.4 are displayed. The clear correlations 
are used to correct the distances and provide a significant reduction in the scat- 
ter. It should be noted that the color has already been corrected for Milky Way 
reddening. 
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Figure 7: Confidence level regions in the (0,m,^a) parameter plane. The 
larger contours show the original results from the High-Z Supernova Survey 
(HZSNS; [15j) and the Supernova Cosmology Project (SCP; [H]). The blue 
areas correspond to the 68.3, 95.4 and 99.7% confidence limits of the Union2 
compilation of SNe la [69]. Baryonic acoustic oscillations (BAO; green) and cos- 
mic microwave background (CMB; orange) constraints from |93] and [Ml are also 
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Figure 8: The current best constraints for a constant w are shown here. A 
cosmological constant {w=0) still is valid solution. The complementarity of the 
different methods is apparent. The supernova probability distribution is nearly 
perpendicular to the ones from the baryonic acoustic oscillations or the microwave 
background. Data taken from [69]. 
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Figure 9: The differences incurred from using different light curve fitters is 
demonstrated here. The same r— and i-filter observations (fluxes, not mag- 
nitudes) of the SNLS supernova 04D2gp (z=0.732) are approximated with the 
SALT2 and the MLCS2k2 fitters, two of the most used fitters at the moment. 
Small deviations in the i— filter fitting can be seen and are due to the different 
underlying templates used in the fitters and differences in how the colors are 
determined simultaneously. 
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Redshift 



Figure 10: Fitted SALT2 color, c, vs redshift for the 3-year SNLS data-set [97]. 
There is strong evidence for a selection effect, with bluer supernovae being found 
at larger redshifts. Since bluer supernovae are intrinsically more luminous this 
effect can be explained as a selection driven by the magnitude limit of the spec- 
troscopic follow-up. 
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Figure 11: Magnification of SNe la in the 3-year SNLS data vs. redshift. Filled 
circles represent the estimated lensing magnification, Am/g^g, to individual Type 
la SNe, based on the best fit halo model from [135) . The effect grows with 
redshift, but the mean value is close to zero, as also indicated by the histogram 
on the side. Figure courtesy of Jakob Jonsson. 




Figure 12: Projected 68% confidence region in the {wo,Wa) plane for a planned 
satelite mission (Eric Linder, private communication) compared with the Union2 
results. While the existing data cannot constrain a time-dependent w future 
missions have the promise to meaningful constrain this quantity. Figure adapted 
by Joel Johansson. 



